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I. INTRODUCTION 

Nonlinear open systems often produce bistabilities or 
dynamical phase transitions. A nice example is the be- 
havior of a Bose-Einstein condensate (BEG) in a pumped 
high-finesse optical cavity, where the nonlinearity is pro- 
duced by the dispersive atom-light interaction [IHS] ■ The 
weak cavity pumping causes a classical electromagnetic 
field to build up between the mirrors. The atoms cou- 
pled dispersively to the radiation field detune the cavity 
according to the overlap of their spatial distribution with 
the mode function of the electric field. Gonsequently, the 
cavity resonance frequency can be shifted away from or 
towards the frequency of the pumping laser; thereby a 
big variation in intensity can be induced merely by the 
spatial redistribution of the atoms. In turn, the inten- 
sity change translates into the variation of the depth of 
the optical dipole potential, and so it acts back upon 
the atomic distribution itself. In a tiny region of the pa- 
rameter settings close to the cavity resonance two stable 
or metastable configurations can exist, giving rise to a 
dynamical phase transition. 

Atom-light interaction itself is a major issue in the 
studies of ultracold quantum gases, being the most uni- 
versal tool in accessing the properties of the system either 
by slowing the cloud of atoms, trapping them in classi- 
cal potentials, putting obstacles in their path, or finally 
detecting and imaging them. Moreover, recent proposals 
have raised the possibility of quantum state preparation 
of the atomic ensemble with the help of measuring the 
output photon signal of a pumped optical resonator (J- 
I13j . The cornerstone of such a quantum state preparation 
with the help of a quantum nondemolition measurement 
is also the mutual backaction between the atomic and 
photonic degrees of freedom. The study of correlations 
between atomic motion and light generated by atom-light 
interaction in an optical cavity is therefore of fundamen- 
tal importance. 

An important research area on the manifestation of 
light-matter interaction is optomechanics, where the ra- 
diation pressure force of a single mode Fabry-Perot res- 



onator is used to manipulate the center of mass motion 
of a mechanical oscillator. For a short review on optome- 
chanical systems and on their experimental realizations 
see Ref. [Mj, and also references therein. The reason 
of popularity of optomechanical studies, besides exper- 
imental realizability, is that the theoretical description 
can be performed relatively simply, involving only the 
few modes of the cavity field and one mode for the mo- 
tion of the mirror |15H17| . Such a paradigmatic system is 
an ideal playground to test correlations between light and 
mesoscopic objects, to understand the underlying physics 
and to speculate on possible applications in quantum in- 
formation processing. 

In recent experiments done with ultracold bosons in 
optical resonators the above concepts unify nicely [THE]. 
The cavity radiation field couples to a single collective 
motional excitation of the Bose condensed atomic sam- 
ple. Starting the experiment with a pure Bose-Einstein 
condensate, other motional excitations can be safely dis- 
regarded and so a situation analogous to optomechan- 
ics can be realized without a movable mirror, but rather 
with the collective motion of an ensemble of atoms. The 
difference between the experimental tools (manipulation 
and detection methods) of traditional optomechanical 
systems and those with ultracold gases nicely comple- 
ment each other, while the theoretical description is very 
similar. 

The aim of the present paper is to discuss correlations 
generated in a hybrid system of ultracold bosons and the 
radiation field, especially close to the region of parame- 
ter settings where the system shows bistability and where 
the optomechanical simplification can be harvested. Pho- 
tons leaking out of the resonator make the cavity field 
noisy, which infiltrates the dynamics of atomic motion 
OUHllTn]. In turn, quantum fiuctuations of the atom 
field have a back-action on the photon statistics. Gor- 
related fluctuations of the light and matter wave fields 
appear then, which are strongly enhanced close to the 
critical regime of bistability. The study of correlations 
is further motivated by the need for justifying the ba- 
sic assumption of the generally used mean-field theories 
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[TBI [^UH^S] . namely, that the atom-photon cross corre- 
lations are negligible and mean values of the light and 
atomic operators can be decoupled. 

The paper is organized as follows. The backbone of 
the paper is Sec. where the model and the theoret- 
ical description of the system in mean-field level is pre- 
sented in great detail. Many of the elements of other 
theoretical models, e.g., the cavity cooling of BEC ex- 
citations [21], spatial self-organization of a BEC in the 
cavity [^, '5^ , and the transient collective atomic recoil 
lasing [27, 28J will be recapitulated here to give a full ac- 
count of the mean field dynamics of a BEC in a cavity. 
The aim, partially, is to reach the optomechanical model 
and discuss the assumptions and approximations to ar- 
rive there. Special attention is paid to the effect of nonlin- 
earities: (i) the nonlinearity caused by atom-light inter- 
action, responsible for the creation of a periodic optical 
potential and also for an effective atom-atom collective 
interaction, and (ii) the nonlinearity caused by atomic 
s-wave collisions. In Sec. |III| we present first mean-field 
results and compare them to experimental observations, 
wherever applicable. Furthermore, the auto- and cross- 
correlations of the quantum fluctuations of the fields will 
be calculated in the stationary state formed by the bal- 
ance of cavity loss and vacuum noise driving. Finally, a 
summary is given in Sec. |IV| 



II. 



THEORETICAL DESCRIPTION OF THE 
SYSTEM 



The system consists of a single mode, high-finesse op- 
tical Fabry-Perot resonator with a waist much smaller 
than the cavity length and a sample of dilute, ultracold 
bosonic atoms prepared to be Bose-Einstein condensed. 
The condensate is supposed to be cigar shaped along the 
cavity axis, with a strong transverse confinement. The 
radiation field inside the cavity is pumped through one of 
its mirrors by a laser with frequency ujp and wavenum- 
ber k = ujp/c, with c being the speed of light. The 
laser frequency is far detuned from the atomic transi- 
tions, so the populations of the electronic excited states 
are negligible. In this limit the atomic internal degrees 
of freedom are frozen; and the atom-light interaction is 
purely dispersive. On the other hand, the cavity fre- 
quency uJc lies close to the pump frequency ujp: the de- 
tuning Ac = i-op — UJc is comparable to k, this latter 
being half of the inverse lifetime of a photon inside the 
cavity. 



A. Hamiltonian 

In the frame co-rotating with the pumping laser field 
the Hamiltonian of the system can be approximated as 



Ha is the Hamilton operator of the ground state atoms 
inside the cavity, given by 

Ha^ f ¥ix)\^^+V,^,{x) + ^¥{x)^{x)\^{x)dx, 
J I 2m ax^ 2 J 

(lb) 

with m being the mass of the atoms, Voxt{x) is the exter- 
nal confining potential along the cavity axis and g being 
the s-wave scattering constant in 1-dimension. The term 
He of the Hamiltonian describes the radiation field of 
the empty, single-mode cavity. 



He = —hAc a^a. 



(Ic) 



The dispersive interaction between the cavity radiation 
field and the atoms in this low excitation limit is given 
by the AC-Stark shift, or light shift: 



H 



hUao^a / ^>\x)'ii{x)cos^{kx)dx, (Id) 



AC 



with Uq being the single atom light-shift, Uq = g^^/A^; 
the unique longitudinal mode function of the single mode 
cavity is cos(fcx) with wavenumber k = LOp/c = 2'k/X. 
The part describing the coupling of the cavity field to 
that of the pump laser is given by 



HcL = V* 



rj a 



(le) 



with 7] being the strength of the driving field; and the 
asterisk stands for complex conjugation. The last part 
of the Hamiltonian, i/vac describes the interaction of the 
cavity field with the broadband reservoir of external ra- 
diation field modes via the partially transmissive mir- 
rors. We will give account for this interaction within the 
Markov approximation, by means of introducing a loss 
rate 2k and a Gaussian white noise operator ^(t) in the 
Heisenberg equation of motion for the field operators [5S] . 



B. Equations of motion 

The equation of motion of the photonic annihilation 
operator is given by 

d r /" 

i—a{t)= -Ac+ •^^{x,t)-^/{x,t)U{x)dx-iK a{t) 

+ iri + i^{t), (2a) 

with U{x) — Uq cos^(fc x) the local single atom light-shift, 
which is a periodic function. Its period is noticeably L — 
A/2, since it contains the mode function squared. The 
operator ^(t) describes Gaussian white noise with zero 
mean and with the only non-vanishing correlation 



{mc\t')) = 2K6{t~t'). 



H = Ha + Hc + Hac + Hcl + H^. 



(la) 



(2b) 



It is nicely exhibited in Eq. ( |2a[ ) , that the dispersive inter- 
action between the atoms and the radiation field causes 
a shift in the resonator frequency proportional to the 
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atomic density "if^ {x , t)'i' {x , t) . However, this frequency 
shift is an operator and couples the equations of motion 
of the radiation field to those of the atomic field opera- 
tors in a nonlinear way. The equation of motion of the 
atomic field operator reads 



iH^^^ix,t) 



2m 



+ Vext(a;) + W(t)a(i)J7(a;) 



+ .g*T(a;,i)$(2;,t) ^{x,t). (2c) 



In the atomic part of the equation of motion (2c 
addition to the inert trap potential V^xt, the atom- light 
interaction creates a A/2 periodic optical potential for 
the atoms, proportional to the dynamical photon number 
operator a^a. 



C. Mean-field solution 

To solve the coupled nonlinear operator equations ([2]) 
simultaneously is a hard task. The most convenient ap- 
proximation is the mean-field approximation, when one 
separates first the operators into a mean value and to 
fluctuations around it: 



*(x,t) 
a{t) 



/Nip{x,t) + S^'{x,t), 
a{t) + 5a{t). 



(3a) 
(3b) 



The mean values are c-numbers, defined as ip{x,t) — 
7V~i/^(vI/(a;, i)) the so-called condensate wavefunction, 
which is normalized to unity; and a(t) — (a(<)) the co- 
herent part of the cavity field. Consequently, the fluctu- 
ations have zero mean. The time evolution of the mean 
values is obtained by substituting ^ into the equations 
of motion ([2| and neglecting all terms containing fluctu- 
ations. By this way one arrives at a Gross-Pitaevskii like 
set of equations of the coupled dynamics 



^-a{t) 



Ac + N{U) 
ffA 



2m 



a{t) + irj, (4a) 
V,^t{.x)+h\a{t)\''U{x) 



-gN\^{xM 



Lp(x,t), 



(4b) 



with the notation ([/) = J tp* {x,t)U{x)(p{x,t)dx. 

There is a simplification of the numerical problem due 
to the possible separation of time scales. The time evolu- 
tion of Eq. ( 4a ) is governed by two characteristic frequen- 
cies, namely the detuning Aq and the photon loss rate k. 



In Eq. ( 4b I the characteristic frequency is set by the recoil 
frequency ujji — hk'^ /{2m). In experiments the latter one 
is usually several orders of magnitude smaller than the 
former ones. For example in the Esslinger group experi- 
ments [IllSlll]? the parameters |Ac| ^ k « 27r x IMHz, 
while the recoil frequency is about loh « 27r x 4 kHz . In 
this situation the dynamics of the resonator field relaxes 
very fast compared to the dynamics of the atomic mo- 
tion and therefore can be considered instantaneous with 



respect to the relaxation of the condensate. One can 
assume for any given atomic configuration that the res- 
onator field has already reached its steady state value, 
which, by Eq. (4a), is 



irj 



N{U) 



(5) 



This steady-state mean field provides the optical poten- 
tial in Eq. ( |4b[ ). So in the end one only integrates solely 
Eq. ( |4b[ ), with a{t) adiabatically eliminated and inserted 
from Eq. ([s]), instead of the coupled Eqs. (4a I and (4b I. 
It is worth noting that if the time scales of the resonator 
field and BEC dynamics do not differ that much, some 
complex coupled solutions can exist [301, which need the 



simultaneous integration of Eqs. (4a) and (4b) 



After the adiabatic elimination of the photon field 
a, the most direct method to numerically calculate the 
steady state of the BEC wavefunction f{x) is the one 
based on the imaginary-time propagation of Eq. ( 4b I . In 



real-time the steady-state solution has the time depen- 
dence: 



Vp(a:,^) = ^(a:)e-'^*/^ 



(6) 



with ^/h the lowest frequency eigenvalue of the nonlinear 
problem \Ah\ . In imaginary time all fluctuations around 
the steady-state die out, since they propagate with higher 
frequencies in real time and consequently vanish faster in 
imaginary time than the steady-state solution. One just 
has to renormalize the solution ^p{x, t) after some finite 
propagation time. Note that, since all quantities on the 
left hand side of Eq. ( 4b ) are real, the condensate wave- 



function (p(x) can also be chosen real. We also note at 
this point that, as we will see later, in the case of an ef- 
fective blue detuning Aq — N{U) > 0, the resonator field 
heats the atomic motion (some excitations have positive 
imaginary parts) and there is no steady-state condensate 
wavefunction at all. However, due to the method of imag- 
inary time propagation, one can find a BEC wavefunction 
even in this case, corresponding to a dynamically unsta- 
ble equilibrium situation. 



D. Fluctuations around the mean-field solution 

Having obtained the steady-state values of the BEC 
wavefunction and the resonator field amplitude, one can 
look for the fluctuations of the annihilation (and cre- 
ation) operators (5a(t), (55'(x,t) {5a) {t), 5'^\x,t)) in lin- 
ear order. This linear stability analysis corresponds to 
the Bogoliubov theory of the BEC system (3lj , and also 
have an analogy in optomechanics, and also in other non- 
linear systems, especially in hydrodynamics. Inserting 
the separation of the field operators ([s]) into Eqs. (2a 



and ( 2c ) , and neglecting fluctuations higher than first 



order, one arrives to 
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i-^Sa(t) = ^ - Ac + N{U) - in 



2™ 



6a{t) + Nass / ip{x)U{x)\d^{x,t) + 6¥ {x,t)\ dx + i^{t) 



S^{x,t) + hU{x)ip{x) a*Ja{t) + assSa\t) 



+ gNip^{x) S¥{x,t) + S'^{x,t) 



(7a) 



(7b) 



r 



where the zeroth order terms cancel, since they fulfill 
Eq. (4b) and ([5| with dass/df = 0. We have introduced 
S^{x, t) = N^'^5^{x, t) e*^*/'^. 

A closer look on Eqs. ([t]) reveals that the time evolu- 
tion of the annihilation operators couple to those of the 
creation operators. It is a consequence of the complex 
nature of the photonic and particle fields. In order to 
diagonalize Eqs. ([t]) one can choose two equivalent ways. 
Either one can separate the complex quantities into real 
and imaginary parts and study their time evolution; this 
is the approach mainly used in optomechanical studies. 
The other way is to diagonalize the set of equations not 
just containing 5a and 5'^ , but also their hermitian ad- 



joints 5a^ and 5¥\ this approach is familiar from the 
Bogoliubov-de Gennes theory of superfluidity. 

Adopting the Bogoliubov-de Gennes way we gather 
the fluctuations into the following column vector R — 
((5a, (5a^, (5^, (5^^)"^, where the superscript T stands for 
transposition, and noises to the other column vector 
Z = (^,^''',0,0)^. Eqs. ([T]) can now be cast into the 
closed form 



ig-^R{t) = MR{t)+iZ{t), 



(8a) 



with 



M 



X ■ fix) 

Y{x) 
A 



A NassX 

-A* -Na*^X 

a;,Y{x) a,,Y{x) fi-^[Ho + gN^^{x)] 

-at,Y{x) -assY{x) -ti-^gNip^{x) 

— + V^^tix) + h\a,,{x)\^U{x) - fi + gN^^{x) 
2m 

Lp(x)U{x)f{x)dx, 

- U{x)(p{x), 

- -Ac + N{U) - iK. 



h-^gNLp^{x) 
h-'[Ho+gNip^ix)] 



(8b) 

(8c) 

(8d) 

(8e) 
(8f) 



Since 6a is not independent of Sa^ and similarly 
is also not independent of d¥ , the matrix M has an 
important symmetry property. It is a consequence that 
the effect of hermitian conjugation of R can be obtained 
with a linear transformation C that swaps the first row 
with the second one and simultaneously the third one 
with the fourth, so i?^^ = C R. It directly follows from 
Eq. (8a) and from this symmetry property that 



CMC 



-M* 



(9) 



In order to study the correlations of the fluctuations 
one has to determine the time evolution of the fluctua- 
tion operators first. One has to introduce quasi-normal 
modes that diagonalize Eq 
simple time evolution. 



and therefore have a 
Let us denote by A'^^ the right 



eigenvectors of M, i.e.. 



(fe) 



(10) 



with Wfe being the corresponding eigenvalue of M. The 
fluctuation operator R can be expanded with the help of 
the eigenvectors (if they form a complete set) 



i?(t)=5]p,(i)rW, 



(11) 



with pk being the operator, or expansion coefficient, of 
the quasi-normal mode k. The operator is given by: 

(12) 

with (•, •) is the Euclidean scalar product, and Z*^*^^ is the 
left eigenvector of M, defined as 



Pkit)^{l^''\R{t)), 



UJ, 



(13) 
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The left and right eigenvectors are normahzed as usual: 
(-;(fe)^^(0) = Sk,i. With the help of Eqs. (pl and 



the normal modes obey the following uncoupled equation 
of motion: 



Pk{t) + iQk{t), 



(14) 



with the transformed noise operator Qk = {l^''\zY On 
integrating Eq. (14), the time dependence of the normal 



mode operators can be obtained: 



-iujkt 



Pfc(0)+ / e"^"'=(*-*''gfe(t'^'. (15) 



Pk{t) 



For a dynamically stable sytem, one needs to have eigen- 
values with negative imaginary parts. In this case fluc- 
tuations (the normal modes) decay to a steady-state. 

The symmetry property ^ has an important conse- 
quence on the spectrum of M. Namely, if a; = e — «7 is 
an eigenvalue of M, than — cj* = —e — i^ is also an eigen- 
value of M. The modes corresponding to these eigenval- 
ues form a pair with positive and negative energies. 

There is an important issue concerning the stability of 
the normal modes: not all of them include the radiation 
components. For simplicity, assume that the external po- 
tential is even with respect to the center of the cavity. In 
this case the combined external and optical potentials is 
also even, and parity is a symmetry of the full system. In 
this case the condensate wavefunction is also even; and 
the matrix M commutes with the parity operator. The 
eigenf unctions of M can therefore be classified by their 
symmetry (being odd or even) . If the condensate fluctua- 
tion parts (third and fourth components) of an eigenvec- 
tor r'^^'^ are odd, than the result of the operator X acting 
on these components of A^'^ is zero, since X contains 
an integration on the whole cavity axis and its kernel is 
even. For these modes there is no coupling term between 
the cavity field and the condensate part [see the first two 
rows of the matrix M in ( 8b |] , and, since damping of the 



modes comes from the cavity decay, these modes remain 
undamped and just marginally stable. Consequently, if 
we want to describe the steady-state values of the correla- 
tions of fluctuations we have to omit those normal modes 
which do not couple to the radiation field (assuming that 
they are initially not populated). 



E. Bloch states, the effects of s-wave scattering and 
that of the collective coupling 



For further analysis we suppose that Vc^ti^) = 0, i.e. 

is the periodic optical po- 



the only potential in Eq. (2c 



tential induced by the resonator field. In this case the 
Bloch functions are good candidates for the complete set 
of single particle wavefunctions: 



(16) 



with Un.q{x) being a periodic function of period A/2. 
Here n is the so-called band index and q G [0,47r/A] is 
the quasi-momcntum of the particle. Af is the constant 
for normalization. If we impose the Born-von Karman 
boundary condition with p periods (the quantization vol- 
ume is pi), then the normalization constant Af ~ P^^^^, 
and the functions Un,q{x) are normalized to unity inside 
a period L = A/2. For practical purposes, we also use 
plane- waves as basis of the Un^q{x) functions (not depend- 
ing on the quasi-momentum q) 



'^n.q {x^ 



1 



(17) 



with ko — 2k = Aw/X. (Of course the basis (17 1 is suit- 
able to express all other L-periodic functions as Fourier 
series.) 

The field operator of the atoms are expanded as 



*(a;) 



'^71, q ' 



(18) 



with br, 



state 



n,q the annihilation operator of the plane-wave 
Here n goes over all integers and q = 
47rm/(Ap), with m £ {0, 1, ... p — 1}. 

The equation of motion of the annihilation operators 
bn,q is obtained easily from Eq. (2c) 



by inserting the field 
operator (|18| and using the orthogonality of the plane- 



wave states. 



^n^^b„^q = 



fj^inko + qY' 
2m 



25 



E 



^nigi ^n2<32^n3i?3'^(n+ni)/co + (g+gi),(n2+n3)fco + (ij2+i?3) ■ (1^) 



The first term corresponds to the kinetic energy, which is 
diagonal for the plane wave states used. The matrix 
appearing in the optical potential is a simple, tridiagonal 



matrix 



(20) 



6 



It is trivial that the kinetic energy is the lowest for 
n = 0, g = 0. The appearance of the optical potential 
makes the Hamiltonian nondiagonal in plane- wave basis, 
however, it mixes only operators with the same quasi- 
momentum q. Coupled plane waves can have indices n 
differing only by An = ±1. On the other hand, the 
s-wave scattering mixes operators with different n and 
different q momenta. The total momentum is conserved 
by the S function imposing {n + ni — n2 — 7^-3)^0 + (? + 
Qi — Q2 — Is) = 0, which allows for normal scattering, 
when q + Qi = q2 + Q3, and n + ni = 71,2 + ^3, and 
also for umklapp scattering processes, where the differ- 
ence of the quasi momentum in the scattering is equal to 
a reciprocal lattice vector: q + qi — q2 + qa — An ■ k^, 
and n -\- ni = n2 + n^ -\- An. It is useful at this point 
to estimate the characteristic frequency corresponding to 
the atom-atom s-wave scattering based on the physical 
parameters relevant to the experimental situation in eg. 
Refs. i, 4J. The s-wave scattering length for the 1) 
states of the ®^Rb atoms is about 5.3 nm. Assuming 
a particle number of = 6 x 10'' atoms distributed 
in the cavity lattice of period L = A/2 = 390 nm and 
the cavity egg-crate potential containing p — 460 val- 
leys, corresponding to a cavity length of 180 /im, and the 
waist of the optical potential as w = 25 /im, the char- 
acteristic frequency of s-wave interaction can be around 
Wsw = 4:TThaN/ (Lpw^m) « 27r x 4 Hz. (Here we have 
neglected any external trapping potential.) This means 
that the characteristic energy of s-wave scattering is 3 
orders of magnitude lower than that of the recoil en- 
ergy separating the bands from each other. As a re- 
sult s-wave scattering dominantly occurs only with atoms 
in the same band and umklapp processes can be ne- 
glected. When incorporating the effects of the confin- 
ing parabolic optical potential in Ref . [5] , the condensate 
density grows roughly by a factor of 100, and results in 
the ratio lOsw/i-^r ~ 0.1. In this situation the contribution 
of umklapp processes is still small enough to expect that 
it does not significantly change the forthcoming results. 

The mean-field equation for the condensate is obtained 
by the substitution 



Vns„ 



(21) 



with /3„ the mean value of the annihilation operators. 
The 6q^o condition for the mean-field comes from the fact 
that the ground state of the translationally invariant sys- 
tem is also invariant under discrete translation, i.e. the 
condensate wavefunction is periodic also with L — A/2. 
(None of the terms in Eq. (21 1 violates the discrete trans- 



q states are not populated for n 7^ 




n = 2 



n = 1 



transitions by s-wave scattering 
= 

transitions by the optical potential 
?i = -1 

n = -2 



FIG. 1: (color online) Schematic diagram of fluctuations cou- 
pling to the the homogeneous BEG in linear order. The grey 
circle represents the macroscopically occupied (BEG) state, 
while the empty circles represent states which are not macro- 
scopically occupied. The arrows show how these states can 
be populated in linear order by scattering from the conden- 
sate via the interaction with the photon field and also via the 
s-wave scattering. 



with ass given by Eq. ([5]), but now with 



([/) =C/o^ /?:(*) W/3„'(i)- 



(23) 



The steady-state components of the condensate ampli- 
tudes f3n{t) evolve as 



/3„(^)=/?„e-^''*/^ 



(24) 



with /i being the chemical potential. It can be seen from 
Eq. (22) that /3n can also be chosen real for all n. Since 
Ln,n' = L-n,-n', and duc to the symmetry of Eq. (22) 



lation of the system, and having q 7^ would cost in 
kinetic energy.) The mean- field equation reads then 



it follows that, /3„ = /3-n, which means, that the Fourier 
expansion of the condensate wavefunction contains only 
cosine terms. 

The equation of motion for the fluctuations, in this 
representation, is given by 



25 

Xp 



n2+n3 J 



(22) 



+ NU, 



[ - Ac 4- N{U) - in da{t) 

+ 5hl^,.)+^m. (25a) 
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25 



n' 



'+"2 I 



ni ,n2 ,n' 



(25b) 



with 



{Ho)n,n' 



ff{nko + qf 



2m 



+ h\ass\'^UoL 

71 



gN_ 
Lp 



'+"2 ~ f^Sn,n'- (25c) 



For a qualitative understanding of Eqs. (251, let us con- 



sider first the weak photon-atom coupling case, when 
either the photon number inside the cavity is small or 
Uq is small. In this situation the atomic distribution 
is essentially unmodified by the periodic potential and 
the condensate wavefunction is almost homogeneous, i.e. 

= Snfil3o. Fluctuations in the atomic density are in- 
duced from the homogeneous condensate. In linear order, 
the atom-light interaction excites fluctuations to bands 
n±l with q — 0, while s-wave interaction populates fluc- 
tuations with arbitrary n and q. However, as discussed 
earlier, in the limit Wsw ^ f^R, s-wave scattering induces 
essentially only intra-band transitions, and therefore, for 
a homogeneous condensate, only states with n = are 
populated by s-wave interaction. Figure [T] depicts graph- 
ically the population of fluctuations excited from the 
homogeneous condensate wavefunction in leading order. 
Moreover, in experiments [3l |4j the one atom light shift 
is a tunable parameter, it's typical value is chosen to be 
on the order of the recoil frequency, i.e. Uq ^ 2tt x A kHz 
and is also much larger than the characteristic frequency 
of s-wave scattering. Therefore on the time scale in which 
the steady-state is reached, the effect of s-wave scattering 
is still not significant. When describing the steady state, 
one can safely neglect it as a first approximation and con- 
sider each band represented by a single state vector with 
q = 0- 

Note that there is an interesting complementary regime 
in which the relevant states are localized, Wannier-type 
ones formed by the coherent superposition of many quasi- 
momentum states with q ^ 0. To describe many body 
effects of the atomic degree of freedom in this regime, one 
can resort to a Bose-Hubbard model with self-consistent 
parameters [5^37| . in which collisions play a vital role 
but intcrband transitions are usually neglected. 



F. Optomechanics 

In the experimental situation of Rcfs. O |4] the BEC 
wavefunction can be considered almost homogeneous, 
with a condensate fraction in the /3o state containing 



6 X 10^ atoms, while the next state with n = ±1, i.e. 
the atomic motional state with wavefunction cos kox con- 
taining only a few hundred atoms. Let us therefore re- 
strict the Hilbert space of the one-particle atomic mo- 
tion into this relevant 2-dimcnsional subspace, containing 
the homogeneous single particle wavefunction and that of 
cosfcoa;. The atomic field operator becomes 



*(a;) ipoix)bQ + ipiix)bi 



(26) 



with the single-particle wavefunctions (po{x) = L 
and fi{x) = y/2jL cos{kQx) , and their corresponding an- 
nihilation operators bo and 61. The single-particle wave- 
function sin(fcox), corresponding to the antisymmetric 
combination of the n = ±1 states, is omitted, since this 
wavefunction doesn't couple to the photonic field due to 
symmetry reasons mentioned above. 

In this model the atomic motion is represented in a 
two-mode Fock space. The mean field expansion is 



(27) 



Nf3,, + 6b„, n = 0,l; 



where /3„ is again the representation of the condensate 
wavefunction in the 2-dimensional Hilbert space, normal- 
ized to unity. The presence of a Bose-Einstein condensate 
distinguishes a subspace of one single mode. Accordingly, 
the fluctuation operators can be expanded to a part par- 
allel to the condensate and to one orthogonal to it: 



Sbn = Pn5b + JnSc, 



(28) 



with 7„ the unit vector orthogonal to the condensate: 
7„ = (— /3o)"^. The part, 5b, parallel to the conden- 
sate can be related to the arbitrariness of the phase of the 
condensate, and corresponds to a zero mode. The orthog- 
onal part, 6c is the sole degree of freedom. Therefore the 
fluctuations of the atomic field, similarly to the photon 
field, has a single degree of freedom. The fluctuations of 
the combined atom-resonator field system has 2 degrees 
of freedom, and is analogous in many ways to cavities 
with a moving mirror, to the so called optomechanical 
systems. 

In order to obtain the mean-field equations and the 
fiuctuation equations in linear order one can start from 
Eq. ([2| and use the truncated field operator ( 26 ) and then 
the mean-field ansatz (27 1. The mean-field equations for 



the condensate now read 



(29) 



with Pi = diag(0, 1), the projection matrix to the sub- 
space of ipi{x); the matrix representing the light-shift is 
given by 



Uk.i = {ipk\Uix)\ipi) 



2 



1 V2/2 
V2/2 1 



(30) 



ass given by Eq. ^, with (U) = Pl{t)Uk,Mt). (Note 
the convention of automatic summation over repeated in- 
dices.) The condensate wavefunction also has the time 
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dependence of Eq. (24), with /3k chosen real. Eq. (29) 
can also be solved either by imaginary time propagation 
or by direct algebraic means using Eq. ( 24 ) and the nor- 
malization condition: /3q + = 1. 

The equations of motion for the fluctuation operators 
can be derived analogously to the way as Eqs. (251 was 
obtained. It is appropriate to separate the trivial time 
dependence due to the chemical potential via the defini- 
tion 6h = N^^^"^ e*''*/'^ 5bk. With this 



^:rM = (-^ - ««)'5a + Nassf3kUk.i{Sbi + Sbj), (31a) 
at 

i^Sbk = Kk.iSbi + {a*Ja + assSa'') Uk,iPi, (31b) 
at 



ith 5^ Ac- N{U), and 



k,l 



'JkU 



-Uk.i 



(31c) 



The matrix K can be thought of being the grand canon- 
ical Hamiltonian of the system, and by virtue of Eqs. 
([29]) and ([24| Kk,i Pi = 0. Equations ( [Ma] ) and ( [31b| ) 
form a linear eigenvalue problem for the fluctuations. It 
is easy to check that the trial function of i5a = and 
Sbk = PkSb is a constant solution, i.e. it is a zero mode. 
It also follows from the normalization of fik that 5b is an- 
tihermitian. With the decomposition 5bk — jS^Sb + 7^ 6c 
one can arrive at a closed set of equations between the 
fluctuations Sa and Sc and their hermitian adjoints. 



i—Sd = {-S - iK)da + Nass{PkUk,ai){Sc + 5c)), 
at 



i-rM = {lkKk,ni) Sc + {■ykUk,il3i) {a*Ja + assSa^) 
at 



(32a) 
(32b) 
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FIG. 2: (color online) The mean cavity photon number 
|qssP as a function of the cavity detuning Ac- The pa- 
rameters are: = 6 x 10'', Uo = 0.96 lor, k — 363.9 cjij. 
The pumping strength is different on the 3 panels: rj = 
(80.06, 283.8, 549.5)aJij for panels a, b, c, respectively. 



The mean-field solution 



Note, that in the optomechanical model s-wave scatter- 
ing can not be included in general, because it populates 
states with g ^ 0, which are disregarded in this model. 
It is a possibility to include some remaining effects of 
the atom-atom collision by coupling only the two states 
under consideration, however, it would result in a nonlo- 
cal and unphysical interaction in coordinate space. Or, 
in order to extend the optomechanical model, although 
still lacking the full consideration of s-wave scattering ex- 
plicitly, the higher quasi-monicntum states can be taken 
into account either by considering a (close to, but non-) 
plane wave condensate wavefunction |25j . or in a two- 
fluid model ESI- 



III. RESULTS 

In the following we summarize and discuss our results 
based on the numerical solutions of Eqs. ^ and ([7|, and 
compare them with those of Eqs. ( [29| ) and ( [32] ) and to 
the findings of the Esslinger group |3j |4j . 



For a comparison between the full GPE solution 
and the optomechanical model we numerically solve the 
Gross-Pitaevskii Equation ( 4b I on a 200 point grid with 



imaginary time propagation with the steady-state value 
of the mean radiation field amplitude ([s]) and compare 
it to the solution of ( 29 1 (also in imaginary time) . The 



results for the mean photon number are plotted in Fig. 
[2] The bistable behavior is nicely exhibited in these fig- 
ures. In panel a) the pumping strength rj — 80.06 
is below a threshold value, the resonance curve gives a 
unique solution for all detunings Ac. Notice, that in 
Eqs. (4a I, or ^ the effective detuning of the resonator 
field is J = Ac - N{U), i.e., the light shift further de- 
tunes the cavity. On the left side of the resonance, the 
effective detuning 5 is negative, the cavity is detuned to 
the red. We will see in the next subsection that this situ- 
ation corresponds to cavity cooling, where the imaginary 
part of the fluctuation spectrum of the combined atom- 
photon system is negative (or nonpositive for the many 
mode system). Fluctuations decay to zero. On the right 
side of the resonance, where 5 is positive, the imaginary 
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part of the fluctuation spectrum changes its sign. Fluc- 
tuations grow exponentially on this side of the resonance, 
that is, the slightly retarded cavity field dynamics heats 
the atomic motion instead of cooling it. Therefore the 
solution above the resonance {S > 0) is dynamically un- 
stable. It is only due to the adiabatic elimination of the 
photon field that the numerical method finds this solu- 
tion as well. Note, that the parameters of our calculations 
mimics those of Ref. @]. The agreement in the photon 
number with Fig. 3 of Ref. [4] is well within the system- 
atic uncertainty of the photon number estimation of the 
experiments (25%). With the use of the recoil frequency 
for Rubidium atoms, Wfl = 2tt x 3.57 kHz, we can com- 
pare the location of the resonance points to the experi- 
mental findings, too. In panel (a), the resonance point is 
at around 28800 ujr w 27rx 103 MHz, in panel (b) the two 
instability points are at around 26700 ujr « 27r x 95 MHz 
and 27700 Wi? w 27r x 99 MHz. All these are within 5% 
of the respective experimental value. In panel (c) the 
resonance points are at around 21500 ojfl « 27r x 77 MHz 
and 27000 cjfl, « 27r x 96 MHz. This latter is only 2% 
higher than in the experiment, but the lower point is 
contrasted to the 27r x 84 MHz value (about 10 % de- 
viation). In our case the bistable region is a bit wider 
than in the experiments. This discrepancy might be at- 
tributed mainly to the uncertainty in the effective de- 
tuning 6 which we obtain by an estimation of the actual 
atom number and overlap between the cavity potential 
and the atomic cloud. 

For ry > r]c{Ac,K,Uo,N), the threshold value of the 
pumping strength depending on the detuning, the cavity 
loss, the light shift and also the atom number, there is 
a region where three solutions exist for the photon num- 
ber and also for the condensate wavefunction. Two of 
them are stable solutions for the numerical method we 
use. These solutions are plotted in panels b) and c) for 
1] = 283.8 ujii and i] — 549.5 ojij, respectively. The un- 
stable solution is not plotted in the figure. We note, 
however, that the upper one of the two plotted solutions 
corresponds to the ^ > 0, cavity heating, situation. So 
this solution is also unstable, dynamically. But since this 
instability is related to photon dynamics (neglected at 
this level) one can still find this solution by integrating 
the GPE for the atoms in imaginary time. 

The fraction of condensate atoms occupying the ho- 
mogeneous and the cos/cq^: states are plotted on Fig. 
[3j The parameter settings are the same as for Fig. [2] 
Both the resonance feature and the bistabily is exhib- 
ited in the condensate component of the cos fcpx mode. 
The condensate is almost homogeneous for the whole 
range, /3j <C /^g ~ 1, and f3f shows a very similar de- 
pendence as the photon number. It is remarkable that a 
very small change in the shape of the condensate wave- 
function causes a drastic increase or fall of the photon 
number. This can be elucidated by the collective cou- 
pling of the atoms to the resonator field. The detuning 
d = Ac — N{U), appearing in the denominator of the 
steady-state value of the resonator field ([5]), can change 
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FIG. 3: (color online) The squares of the condensate com- 
ponents /3o and /3f as a functions of the detuning Ac in a 
semi-log scale. The parameter settings are the same as for 
Fig-i 



a lot even if the variation of (U) is small, because of the 
large number N multiplying it. 



B. Fluctuation spectrum 

The linear stability analysis of the mean-field solution 
is done with the help of the linear equations ([T]) for the 
case of the full model, and with the help of Eqs. (32) 
for the optomechanical approximation. Once the con- 
densate wavefunction and the steady-state value of the 
mean cavity field is known, one can construct the ma- 
trix M and calculate its eigenvalues and eigenvectors. In 
the full model, discretized on a grid with 200 sites, M 
is a 402 X 402 matrix, since it acts on a row vector hav- 
ing 2 elements for the photon fiuctuation operator and 
its hermitian adjoint, and also having 400 elements for 
the discretized field fluctuation operator and on its her- 
mitian adjoint. For the optomechanical model one has 
a much smaller, 4x4 matrix, since the zero mode is al- 
ready separated in this case, so one has 2 components for 
the photon fluctuations and 2 components for the atomic 
motion orthogonal to the condensate wavefunction. The 
numerical diagonalization was obtained with the help of 
the LAPACK package. 
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FIG. 4; (color online) The fluctuation spectrum oj^ — + i'^k of the atom-cavity dynamics. The real and imaginary parts of 
the eigenvalues are plotted as functions of Ac- AH quantities are measured in units of the recoil frequency. The parameter 
settings are the same as for Fig. [2] Note, that the fluctuation spectrum of the full model, discretized on a 200 grid, contains 
201 pairs of eigenvalues. We only plot those 2, which correspond to those of the optomechanical system. 



Note that M is a general complex matrix. Its eigenval- 
ues are complex. However, due to the symmetry relation 
([9]) the eigenvalues of M come in pairs. Each element 
of the pair have the same imaginary part, while the real 
parts are just the opposite of each other. Figure |4] shows 
the eigenvalues of M for the same parameter settings as 
that of Figs. [2] and [3] The left panels show the real part 
Ck of the /cth eigenvalue, while the right panels show the 
modulus of the imaginary parts 7^ for uj}~ = e^+i^k. Only 
those two eigenvalues are plotted, which are present in 
the optomechanical model and have positive real parts. 
The other two eigenvalues of the optomechanical model 
can be obtained simply by changing the sign of the real 
parts to negative: Wfc — >■ — efc + i^k- In the bistability 
regime, the excitation frequencies are presented for both 
mean field solutions. 

Again, all figures show that the optomechanical ap- 
proximation very well reproduces the results of the full 
GPE simulation even for dynamical quantities. The 
two modes are easy to be physically interpreted. The 
first eigenvalue, having a bigger real part around 5 = 
Ac — N (U) and an imaginary part almost exactly at — k, 
corresponds to a mainly photonic mode. Note that the 
imaginary part of this mode remains also negative above 



the resonance. The other eigenvalue, having a real part 
around -iujfi corresponds mainly to atomic fluctuations. 
The imaginary part of this eigenvalue changes sign at res- 
onance. Below resonance the imaginary part is negative 
(fluctuations in the atomic motion are damped due to the 
interaction with the resonator field [331110]), while above 
resonance the imaginary part becomes positive (fluctua- 
tions in the atomic motion are exponentially growing in 
time). The situation is clear on panel a), where the reso- 
nance point can be exactly defined. On panels b) and c) 
the coexistence of the two solutions makes the definition 
of the resonance point ambiguous. Nevertheless, on pan- 
els b) and c) the curves which continue to the left hand 
side of the figure correspond to the solutions with cavity 
cooling and those which continue to the right hand side 
correspond to the heating solution. The dynamical cool- 
ing and heating effects are closely related to the same 
effects at the single atom level |1T| - H5] . 

In Refs. [21 H] the coherent atomic dynamics was 
also studied. The harmonic oscillator behavior of the 
low-energy atomic dynamics caused periodic oscillations 
in the output photon signal with a frequency close to 
35 kHz, as reported and explained in Ref. [3], or close 
to around 42 kHz in Ref. [4j. In our model the dynam- 
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ics of the two coupled harmonic oscillators of the res- 
onator field and the atomic collective motion is described 
by Eqs. (32). These are linearized equations of motion 
around the steady-state configurations. For very small 
fiuctuations the semiclassical time evolution can be in- 
terpreted as orbitals around the steady-state fixed points 
with frequencies plotted in the left panels of Fig. [4] Due 
to the imaginary parts the trajectories spiral closer to 
(for fluctuations around the cooling solution) or farther 
from (for fluctuations around the heating solution) the 
corresponding fixed point. From Fig. [4] we can see, that 
the frequency of atomic motion dominated fluctuations is 
close to 4 LUn for parameters not very close to the bistable 
region. In the bistable regime two steady-state solutions 
exist, and correspondingly there are two flxed points in 
the harmonic oscillator phase space. According to Fig. 
|4]the oscillation frequency for the cooling fixed point re- 
mains around 4aj/j and softens close to the endpoint, 
while for the heating solution the oscillation frequency 
grows to around lOw/j before softening at its endpoint. 
This kind of renormalization of the atomic dynamics by 
the interaction with radiation is often referred to as the 
'optical spring effect'. The approximately 10 wj^ angular 
frequency quantitatively gives back the experimentally 
found 42 kHz oscillation frequency of the Fig. 3C inset 
ofRef. g. 

It is interesting to point out that in the bistable region 
even small atomic fluctuations can result in a phase-space 
trajectory that orbits around both of the fixed points 
[see Fig. 3 of Ref. [3]]. In this case such a lineariza- 
tion strategy simply can not work because the transition 
between the two fixed points is necessarily a nonlinear 
effect. However, if the trajectory would lie in the basin 
of the linear region of the corresponding fixed points ev- 
erywhere except the small region of the separatrix, one 
could expect the oscillation to be described by the above 
two frequencies: when the system is in a part of the tra- 
jectory inside the attraction basin of the cooling fixed 
point the angular frequency would be around 4a;i{, while 
after crossing the separatrix it would change to around 
lOojfl,. Such a trajectory of coherent oscillations would 
also cause a periodic output photon signal but both of 
the frequencies would appear, i.e. the count rate peaks 
would come in pairs. 

Figure [S] depicts the regions in the parameter space 
which represent the heating and (cooling) solutions, i.e. 
the solutions with dynamical instability (stability). The 
bistable region, where both a heating and a cooling so- 
lution can be found, is wedged between the two distinct 
regions. The tip of the wedge corresponds to the critical 
point {AccVc)- 



C. Correlations and entanglement 

Second order correlations of the fluctuations are calcu- 
lated with the help of the quasi-normal modes and the 
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FIG. 5: (color online) The nonequilibrium phase diagram of 
the system indicating the parameter regions with full dynam- 
ical stability (the cooling region) and the region with a dy- 
namical instability caused by the energy transfer of the cavity 
(the heating region). The horizontal lines represent the pa- 
rameter r] at values corresponding to the subfigures a) b) and 
c) of other figures with numeric results. 



noise correlations, since by virtue of Eqs. (11) and (151 
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(33) 

For a stable system, where all the eigenvalues have nega- 
tive imaginary parts, the first term on the right hand side, 
corresponding to the initial condition of the fluctuations, 
vanishes for times much longer than the characteristic 
decay times of the system. For the second term one can 
use the noise correlation function (2b), the definition of 
the vector noise Z and that of Q. For t — >■ c» it follows 
straightforwardly, that 



{Rk{t) Ri{t)) 



;(»»)* ;(")* „("') 

2^E V J ? ' (34) 



where we index the components of the row vectors R and 
Z starting from 1, and also dropped the exponential term 
vanishing for large times in the case of a stable system. 
Note that the zero mode, i.e., Sb in (28) representing the 



phase fluctuations of the condensate, do not contribute to 
this sum because its eigenfunction does not have photon 



component, i.e., ll 



(zero mode) i(zcro mode) 

— '9 



= 0. 



To be able to relate our results more explicitly to other 
works, let us introduce quadrature operators, according 
to 5x = {5a + 6a^)/V2, 6y = --i{Sa - Sa^)/V2, 6X = 
{Sc + Sc^)/V2, 5Y = -i{Sc-6c'')/V2. These quadratude 
operators are hermitian operators, easily expressed with 
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FIG. 6: (color online) The nonclassical part of the photon 
number vs. the cavity detuning of the optomechanical model. 
For reference we have also plotted the mean-field solution, 
|qss|^, with a dashed line. All parameters are the same as for 
Fig-i 



FIG. 7: (color online) The depletion of the condensate vs. the 
cavity detuning of the optomechanical model. All parameters 
are the same as for Fig. [2] 



the help of the field operators R. We assemble them into 
the following row vector: u — {Sx,Sy,6X,SY)'^ . With 
the quadratures being hcrmitian, one can define a real 
correlation matrix by 

CkAt) = ^{uk{t)ui{t)+ui{t)uk{t)), (35) 



which is in the following block form 
C 



P X 
A 



(36) 



with P representing the correlations of the photonic de- 
gree of freedom, A representing atomic fluctuations, and 
X describing the cross correlations. For example, on top 
of the mean field contribution |ap, there is a nonclassical 
part of the photon number which is given by 



2,2 



1 



n' 1-1 = {So' So) = 

(37) 

In an empty resonator close to zero temperature, where 
{£,Ht)£^{t')) — 0, the nonclassical part of the photon num- 
ber (|37|) vanishes, and the resonator field is in a pure 



coherent state. The fluctuation of the quadratures are 
distributed equally and the Heisenberg uncertainty prin- 
ciple is fulfllled in a sharp sense ((<5x^) = {Sy"^) = 1/2). 
However, due to atom-photon interaction, the photon 
fluctuations couple to those of the atomic motion. By 
iterative substitution of Eqs. (32) one can simply check 



that Tipij = (Sa^Sa) is no longer zero in the presence of a 
Bose-Einstein condensate. The photon field is no longer 
in a purely coherent state, and the correlation matrix P 
is not simply half of the unit matrix. Indeed, since Sa 
couples now to Sa^ (through Sc) and the coupling term is 
proportional with a^g, which is complex, the correlation 
matrix P is not isotropic. The angle of the major axis 
coincides with twice the phase of ass, taking the values 
from — TT (when the system is far from the resonance) 
to zero (at resonance). By performing a rotation of the 
correlation matrix by this angle to bring it into a diag- 
onal form, one of the eigenvalues will remain 1/2, as in 
the pure coherent-state case, while the other one will al- 
ways become bigger than 1/2. Since the trace of a matrix 
is invariant under rotations, the nonclassical part of the 
photon number n'^^^ = (A> — l/2)/2, with A> being the 
bigger eigenvalue of P. Figure [6] shows the nonclassical 
contribution to the photon number. The classical average 
I ass I of Fig. [2] is also plotted here with a dashed line for 
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FIG. 8: (color online) The logarithmic negativity of the con- 
densate vs. the cavity detuning of the optomechanical model. 
All parameters are the same as for Fig. [2] 



depletion scales completely differently than that caused 
by s-wave scattering and is strongly influenced by the 
presence of the resonance. The diffusion of atoms out 
of the condensate can also be interpreted as a quantum 
measurement-induced back action process which stems 
from the dispersive atom-light interaction [5j I19j and oc- 
curs also even in phase contrast imaging of a conden- 
sate where the photon field is propagating in free space 
[ini HZ] ■ In a recent paper [TH] we have shown that the 
depletion has a large steady-state value even in the limit 
of vanishing interaction strength Uq. For a fixed large 
detuning, Ac 3> k, NUq, the amount of non-condensed 
atoms was estimated by Ac/wr, and is connected to the 
ratio of the photon energy and the energy of motional 
excitations. In Fig. [7] the detuning is a variable and the 
above condition is not fulfilled. However, far from reso- 
nance, the same ratio determines the depletion with the 
effective photon energy given by the detuning 6c- 

The amount of entanglement between the atomic mo- 
tional and photonic degrees of freedom can also be calcu- 
lated with the help of the correlation function assuming 
that the state of the system is a Gaussian one. The log- 
arithmic negativity, E_\f is a useful measure of entangle- 
ment in our case, since it can be directly calculated with 



the help of the correlation matrix ( 35 1 



E_\f = max(0, — In 2r/ 



(39a) 



where 



2-^/^Jj:{C) - v/S(C)2 -4detC 



(39b) 



reference. It can be seen, by comparing the two curves, 
that the nonclassical contribution to the photon number 
is very small compared to the contribution of the mean 
field far from the resonance. When approaching the res- 
onance, the nonclassical part grows faster than the mean 
field part, and, in the vicinity of the resonance, it even 
exceeds the mean field contribution. The nonclassical 
contribution of the photon number is plotted only for 
the cavity cooling regime, where a steady state solution 
exists. 

Similarly, the number of atoms outside the condensate, 
i.e., the depletion, is evaluated as 



N' = ((5c^ 5c) = NiSc^ Sc) = 



C4 , 



1 



(38) 



Figure [7] shows the steady-state number of particles out- 
side the condensate for the parameter settings of the ear- 
lier plots. Notice that on the right hand side of the reso- 
nance depletion is not plotted. In this regime, where the 
cavity heats atomic motion instead of cooling it, there is 
no steady-state condensate, and depletion grows in time 
exponentially. The steady-state depletion is analogous 
to the quantum depletion of the ground state of a non- 
ideal Bose gas of atoms due to collisions. In the cavity 
case the interaction between the atoms is provided by 
the collective coupling to the photon field. However, this 



is the smaller symplectic eigenvalue of the two-mode 
Gaussian state, with S(C) = det P -|- det A — 2detX. 
The state is an entangled state if and only if Ej^ ^ 0. 
The larger the logarithmic negativity, the larger the en- 
tanglement between the atomic motion and the photonic 
degree. 

Figure |8] shows the steady-state value of the logarith- 
mic negativity, Ej^ as a function of Ac for the three 
parameter settings of Fig. [2] These results hold only for 
the cooling solution. In the heating regime, where there is 
no steady-state, the entanglement between the photonic 
mode and the atomic motion also grows in time. The log- 
arithmic negativity takes very small values in the whole 
range of the presented parameters except for a very small 
region around the instability. Apart from this narrow re- 
gion, the entanglement is small even compared to the val- 
ues of other optomechanical systems 17! . The smallness 
of the entanglement might be attributed to the big dif- 
ference between the occupation numbers of the photonic 
and atomic modes, or equivalently, to the large difference 
in the effective energies of the decoupled subsystems. To 
reach higher values of entanglement either the timescales 
of the photonic and atomic degrees of freedom should be 
closer to each other, or a much stronger driving is needed 
to attain high photon numbers. 
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IV. DISCUSSION 

In this paper we have investigated the one-dimensional 
dynamics of a Bose-Einstein condensate inside a driven 
optical cavity. As the dispersive atom-photon interac- 
tion couples the atomic motion to the dynamics of the 
photonic field in a nonlinear way, strong correlations can 
appear. The strength of this coupling is inversely pro- 
portional to the detuning of the pump frequency from 
the atomic transition, therefore it can be tuned in ex- 
perimental implementations. We recited and analyzed 
the mapping of the original system to a two mode ef- 
fective model in which only the two highest populated 
one-particle states are kept from a plane-wave expan- 
sion of the atomic motion [SHS] . By solving the coupled 
Gross-Pitaevskii equations it was possible to reproduce 
the bistable behavior caused by the nonlinear coupling 
[U m and to provide a phase diagram of the system 
partitioning the whole parameter space into regions with 
full dynamical stability of the mean-field solution, the 
cooling region; to a region with a dynamically instabil- 
ity attributed to cavity heating; and to a region where 
both a stable and an unstable solutions can exist. We 
have compared the mean-field solution and the fluctu- 
ation spectrum of the optomechanical model to that of 
the model not restricted to the first two highly occu- 
pied modes. In the cavity heating region, the unstable 
polariton mode can have a positive imaginary part in 
the order of a kilohertz giving an evaporation rate of the 
Bose-Einstein condensate in milliseconds. Such timescale 
is in the experimental reach. 

The dispersive atom-photon interaction not just causes 
cavity cooling or cavity heating but also alters the statis- 
tics of the constituent subsystems. In the framework 
of the optomechanical model, second order correlations 
were also investigated between the radiation field of the 
cavity and the motional mode of the Bose-Einstein con- 



densate in the cavity cooling regime. Significant contri- 
butions beyond that of the mean field were found to the 
photon and particle numbers. The strong depletion of 
the condensate shows some analogy with the excess noise 
in lasers ,48. ,49j , already discussed in our previous work 
|18j . It is interesting however, that the huge nonclassical 
contribution in the autocorrelation of the photonic and 
atomic operators do not manifest in the entanglement 
of these variables. The lack of entanglement can be at- 
tributed to the big difference of the occupation of these 
modes. While the atom number was assumed to be in 
the order of 10^, the photon number ranged in the order 
of unity. 

The experimental progress in combining cavity QED 
systems with ultracold atoms promises an interesting 
playground to test the manifestation of light-matter in- 
teractions in the mesoscopic scale. In such systems both 
the radiation and the atomic part are dynamical enti- 
ties. The better understanding of their interplay can 
have impact not just on our knowledge of nonequilib- 
rium systems, but also on implementations of quantum 
information processing devices or quantum simulators of 
other systems. 
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